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Merger and non-merger galaxy clusters in 
cosmological AMR simulations. 
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Abstract. Aims: We discuss the dependence of shocks, cosmic rays acceleration and turbu- 
lence on the dynamical state of the host clusters. Method: We perform cosmological sim- 
ulations with the grid code ENZO 1.5, with a mesh refinement scheme tailored to follow 
at high resolution shocks and turbulence developed in the clusters volume. Results: Sizable 
differences are found when some important properties, connected to non-thermal activity in 
clusters, are compared for post-merger, merging and relaxing systems. 
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1. Introduction 

Numerical simulations presently provide a 
unique way to study the generation and evo- 
lution of shock waves, turbulence and chaotic 
motions following the evolution of large scale 
cosmic structure over a wide rang e of scales, 
and in a fully-time dependent way (Dol ag et alJ 
2008 and references therein). Simulations can 
also be a very powerful approach to tackle the 
problem of explaining the occurrence and mor- 
phologies of non-thermal emissions in clusters, 
obser ved in many systems (e.g. iFerrari et alJ 
|2008| for a review). We recently employed a 
tailored Adaptive Mesh Refinement (AMR) 
scheme in the ENZO cosmological code ( e.g. 
lO'Shea etalJ 12004 iNorman etalJ I2007I) in 
order to study these phenomen a with un- 
precedented high dynamic range (IVazza et alJ 
l2010ah IVazzaetal.ll2010dl) . Here we summa- 
rize some of the most important findings of our 
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Fig. 1. Map of gas temperature for a major 
merger at z ~ 0.6. The side of the image is 
8.8 Mpc/h, the depth along the line of sight is 
25kpc/h. 
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work, reporting dependences between shocks, 
CR acceleration and turbulence and the dy- 
namical state of the simulated host clusters. 

2. ENZO Simulations. 

The computations presented in this work were 
performed using the ENZO 1.5 code devel- 
oped by the Laboratory for Computational 
Astrophysics at the University of California 
in San Diego (http ://lc a.ucsd.edu), s ee also 
lO'Shea etall d2004l) and lNorman et all d2007l) . 

We performed non-radiative ACDM sim- 
ulations sampling a total cosmological volume 
with the size of Lbox ~ 440 Mpc/h, in which we 
re-simulated the evolution of the 20 most mas- 
sive galaxy clusters, with a DM mass resolu- 
tion of ;«dm = 6.76 • 10 8 Mq and a peak spatial 
resolution of Ax » 25 kpc/h in side < 3R V u- for 
each cluster dVazza et al.l2010al) . Our Adaptive 
Mesh Refinement (AMR) strategy was tailored 
to keep the maximum available resolution also 
at large distances from the cluster centers, 
tracking the propagation of strong discontinu- 
ities in the velocity field ass ociated with shocks 
or turbulent motions (see IVazza et al.1 l2009bl 
for a detailed discussion). At z — each sim- 
ulated clusters is sampled with a high number 
of cells, Ngrid ~ 500 3 - 600 3 , allowing a study 
of the intra cluster medium (ICM) across 2-3 
orders of magnitude in spatial scales (see Fig. 
[T]for a representative example). 

Our clusters have total masses in the range 
6 • 10 14 < M/Mq < 3 • 10 15 ; they were di- 
vided according to their dynamical state, fol- 
lowing in detai l their matter accre tion history 
forz < 1.0 (see lVazza et al.ll2010al for details). 
According to this definition, our sample con- 
tains 10 post-merger systems (i.e. clusters with 
a merger with a mass ratio larger than 1/3 for 
z < 1), 6 merging clusters and 4 relaxed clus- 
ters at z = 0. 

3. RESULTS 

3.1. Shocks and Cosmic Rays 

We identified shocks in t he ICM with t he sam e 
procedure presented in IVazza et ail (|2009a), 
based on the analysis of velocity jumps across 




Fig. 2. 2-dimensional slice showing the Mach 
number of shocks for cluster El of our sample. 
The side of the slice is 1 3 Mpc/h and the depth 
along the line of sight is 25 kpc/h. 



close cells (see Figure [2] for a map of the re- 
constructed Mach numbers). All clusters show 
similar shocks frequencies and energy dis- 
tribution inside R v [ r : the profiles of average 
Mach numbers are very fiat and character- 
ized by very weak shocks, M < 2 (top 
panel of FigfSJ. Sizable and stronger shocks 
(M ~ 2.5 - 3.5) are detected only in one 
post-merger systems and in one merging sys- 
tem. The volume distribution of shocks is very 
steep (a < -4 with a = d log N( M)/d log M), 
while the distribution of flux of thermal energy 
through shock surfaces shows a well defined 
peak of thermalization at M ~ 2 (see bottom 
panel of Fig [3]). In general our distributions are 
slightly steeper than what w as previously re- 
ported in the li t eratur e (e.g. iRvu et al.l l2003t 
iPfrommer et"ai1 120061: ISkillman et alj boOSl) 
which can be explained as an effect of differ- 
ent resolution, re-ionization models and shocks 
detection schemes. Only for the few M > 2 
cases reported above the flux distribution is 
flatter, with a ~ -2. We investigated the in- 
jection of relativistic protons (CR) at shocks 
applying the e f ficien cy function introduced by 
iKang & Jon es (2007) in the framework of the 
Diffusive Shock Acceleration model. We found 
that the ratio between the CR energy flux and 
the thermal flux is < 0.05 inside Q.2R v [ r and 
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Fig. 3. Top: radial profile of average Mach 
number for the shocks in our sample. Bottom: 
distribution of thermal energy flux at shocks in 
the cluster samples. 



< 0.1 inside R vlr . Only in one strong post- 
merger system we measure a CR injection of 
~ 0.2 — 0.3 of the thermal energy flux. The 
estimated amount of CR energy inside clus- 
ters implied by our results is well below ~ 
0.1, which is presently allowed by observa- 
tional upper limits coming from the non detec- 
tion of gamma radia tion from secondaries (e.g. 
lAleksicetalJ 120101) and from the non detec- 
tion of diffuse radio emi ssion in relaxed clus- 
ters dBrunetti et al.ll2007l) . 



3.2. Turbulence 

The turbulence in the ICM is expected to be 
sustained by the hierarchical process of mat- 




Fig. 4. Radial profile of the ratio between tur- 
bulent energy inside a given radius, and the 
thermal energy. The colors are as in Fig [3] 



ter accretion in evolving galaxy cluste r s (e.g 



Norman & Brvanl Il999t iDolag et all 120051: 



Iapichino & Niemever! 120081) . In Vazza et al.1 



(2010d) we discussed the application of 2 dif- 
ferent filtering techniques to disentangle the lo- 
cal mean velocity field of the gas (assuming it 
is laminar above a certain scale) and the tur- 
bulent velocity field. Despite small differences 
from case to case, the two methods provided 
consistent estimates of the turbul ent energies 
of our clusters (see Sec. 3.1 of IVazza et aU 
2010d). In Fig0] we show the radial profile 
of the turbulent to thermal energy ratio (for 
the filtering at the scale of « 300 kpc) for 
all clusters in our sample. As expected turbu- 
lent motions are strong in merging and in post- 
merger systems (-EWb /^therm ~ 0.1 — 0.3 for 
r < 0.2R wlr ), while they are rather weak within 
the cores of relaxed systems CEturb/^therm < 
0.05). However, almost all clusters host a siz- 
able amount of turbulence (~ 20 - 40-Ethemi) 
within the total virial volume. An interesting 
finding of our simulations is that, on average, 
the turbulent energy present in merging sys- 
tems is slightly larger than the turbulent en- 
ergy of post-merger systems, when normalized 
to the thermal energy of the host cluster. This 
happens because the thermal energy inside the 
two clusters is initially lowered as an effect 
of the accretion of cold gas from filaments, 
while the turbulent energy almost stadely in- 
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creases (see panels of Fig[5] and also th e dis- 
cussion in Sec. 3.2 of lVazzaet a fll2010d). Our 
sample is large enough to allow the statisti- 
cal study of the occurrence of turbulence in 
cluster cores, and to compare with some basic 
expectation from th e theoretical "turbu lent re- 
acceleration" (e.g. iBrunetti et al]l200ll) . as for 
instance with the estimated frequency of "tur- 
bulent" clusters. We investigated this issue in 
detail, by considering the typical volumes and 
turbulent energies needed to produce realistic 
radio halo emission: at z ~ 0, about 1/3 of 
our simulated clusters host a level of turbu- 
lent energy of ~ 0.25£' t herm inside a volume 
of ~ /?vir/3, in line with existing observations 
(Fig|6]l. Despite these rather large values of tur- 
bulence in merging or post-merger clusters, the 
ICM at the smallest scale is only weakly turbu- 
lent, as follows from the general shape of the 
velocity power spec t rum of our simu lated ICM 
dVazza ei all 1201 Oct IXu et al.ll2009h . If turbu- 
lence in our clusters is computed for the typical 
spatial scales of real X-ray obs ervations per- 
forme d with XMM-Newton (e.g. lSanders et al.l 
2010), ~ 30 kpc, the simulated turbulent ve- 
locities of the ICM within the core region of 
clusters are much smaller than the upper limits 
from observations (see FigJTJ. 

3.3. Azimuthal Scatter of Gas Profiles 

We investigated the level of intrinsic azimuthal 
scatter from the center of our clusters, fo- 
cusing on some important thermodynamical 
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Fig. 6. Cumulative distribution function for 
^tui-b/^therm inside R vir /2, R vil /3 and R vil -/IQ 
for the simulated clusters assuming I < 
3Q0kpc for turbulence. The solid line shows 
the differential distributions for R vu /3. The 
grey band shows the turbulence required in 
the turbulent re-acceler ation scenario (e.g. 
IBrunetti & Laz arian 201(J. 



quantities (e.g. density, temperature, entropy 
and X-ray luminosity). This is interesting be- 
cause recent long SUZAKU exposures mea- 
sured sizable differences between the clus- 
ter pro files in different sectors from the cen - 
ter (e.g. lBautz et"aHl2009tlGeorge et al.ll2009h . 
The inspection of the radial properties of our 
simulated data offers a way to assess the need 
of non-gravitational physics in simulations, 
such as AGN, Cosmic Rays, magnetic field, in 
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T[keV] 

Fig. 7. Relation between temperature and tur- 
bulent velocity dispersion (small squares) for 
our clusters (squares in colors ) and for the 
XMM -Newton observations of ISanders et alJ 
(2010). The thick squares are for turbulent field 
below < 30kpc. 



order to reconcile with observations. In Vazza 
et al. (2010 c) we tackled this issue by analyz- 
ing ENZO and GADGET2 run s of comparable 
resolution (Dol ag et al.l 12005). We processed 
all cluste rs with the filter i ng tec hnique intro- 
duced in iRoncarelli et al.l d2006l) . designed to 
remove the densest clumpy gas component 
from each cluster dataset, and computed the 
cluster profiles in sectors of different angular 
size. In Figj8]we report the azimuthal scatter 
(Aar = 7r/8) for the different quantities in the 
ENZO clusters. Post-merger clusters present 
the largest degree of azimuthal scatter at all 
radii while the relaxed ones have the smallest 
scatter. The scatter increases of a factor ~ 2 
going from R vlI to 2 • R vn ; the scatter in tem- 
perature is always larger than the scatter of the 
other quantities. Once that the contribution of 
dumpiness is removed from the datasets, post- 
merger clusters are the one more affected by 
an asymmetric distribution of the thermal gas, 
due to the chaotic pattern induced along merg- 
ers and filamentary accretions. The amount of 
azimuthal scatter is fou nd to be in fine with ex- 
isting observations (see IVazza et al.ll201 0b for 
a detailed discussion), even if the presence of a 
un-removed clumpy component in the observa- 
tion should be assumed to fit the case of clus- 




Fig. 8. Azimuthal scatter profile for our sam- 
ple, after the "99 percent" filtering of the ph- 
syical fields. The different colors refer to the 
different physical quantities, the linestyles re- 
fer to the dynamical classes of the sample. 



ter PKS0745-191 dGeoree et alJl2009t> . These 
results suggest that non-gravitational processes 
may be not needed in simulations to model the 
thermodynamics of t he gas in the outer clus- 
ter re gions (see also iBurns et alj|20l"ot iNagall 
120111) . 



4. Discussion and Conclusion. 

In summary, our AMR simulations with ENZO 
show that, for all (or most) of clusters: 

- the volume and energy flux distribution of 
shocks in the ICM are steep and most of 
shocks have M < 2; 

- the radial distribution of Mach numbers is 
very flat inside R v i T ; 

- the average CR injection efficiency injec- 
tion is small, ecR/ethenn < 0.05; 

- the kinetic energy of turbulent motions is ~ 
0.2 - 0.4 of the total thermal energy inside 
^vir; 

- the turbulent motions of the ICM at scales 
< 30 kpc are subsonic v tur b < 0.1c c ; 

- the azimuthal scatter of the clusters profiles 
increases by ~ 2 from R wn to 2R WU . 

We find, however, the following significant 
dependences on the dynamical state of clusters: 
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- strong M ~ 2.5-3.5 shocks are found 
only inside the central region of clusters in 
~ 1 / 10 of cases (never in relaxed clusters); 

- in these rare systems, the injection effi- 
ciency of CR can be quite large, ecR ~ 

0.1 - 0.2etherm; 

- in post-merger system, the turbulent energy 
inside R v i T /3 is E tul \, ~ 0.25E± erm , while it 
is ~ O.lZstherm in post-merger systems and 
~ 0.05£therm in relaxed ones; 

- post-merger systems present a larger vol- 
ume filling factor of turbulent motions 
compared to relaxed systems, and host 
enough turbulent energy to explain radio- 
halos within the re-acceleration scenario 
(see Brunetti et al. 2007); 

- the azimuthal scatter in merging and post- 
merger clusters is larger (by a 20 - 50 per- 
cent) compared to the relaxed ones. 
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